Phase Diagram of the Attractive Hubbard Model with Inhomogeneous Interactions 
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The phase diagram of the attractive Hubbard model with spatially inhomogeneous interactions is 
obtained using a single site dynamical mean field theory like approach. The model is characterized 
by three parameters: the interaction strength, the active fraction (fraction of sites with the attrac- 
tive interaction), and electron filling. The calculations indicate that in a parameter regime with 
intermediate values of interaction strength (compared to the bare bandwidth of the electrons), and 
intermediate values of the active fraction, "non-BCS" superconductivity is obtained. The results of 
this work are likely to be relevant to many systems with spatially inhomogeneous superconductivity 
such as strongly correlated oxides, systems with negative U centers, and, in future, cold atom optical 
lattices. 
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I. INTRODUCTION 



There is a growing body of experimental evidence that 
strongly correlated oxides such as cuprateSfi'^'^'^'^'^'i^i^ 
manganite a^°i^^'^^i^^'^'^i^^ etc., are electronically inhomo- 
geneous. The term "electronically inhomogeneous" is 
used to describe states with spatially inhomogeneous 
electronic orders, i.e., of orders of charge, spin, super- 
conducting gap etc. There are suggestions that such in- 
homogeneous electronic states are one of the character- 
istic features intrinsic to strongly correlated materialsi^ 
arising out of their "electronic softness" A clear un- 
derstanding of this phenomenon could, for example, sug- 
gest possibilities of controlling the nature and size of the 
electronic inhomogeneities, and can lead to, inter alia, 
possible device applications of these materials. Efforts 
directed towards uncovering the physics of the origin and 
nature of electronically inhomogeneous electronic states, 
therefore, have emerged as a very active research area. 

Of particular interest to this work is electronically 
inhomogeneous superconducting state which is found 
in many systems of current interest. High temper- 
ature superconducting cuprates arc prominent exam- 
ples of systems showing inhomogeneous superconduc- 
tivity. The past five or so years have witnessed fas- 
cinating experimental work based on scanning probes 
that have revealed a wealth of information regarding the 
nature of the inhomogeneous superconducting state in 
cuprates i^i^ii^i^i^iii^i^ In particular the experimental work 
reported in references has clearly demonstrated a 
distribution of gaps, and even regions with gaps above the 
superconducting transition temperatures. There could 
be several physical origins to this phonemenon, such 
as one body disorder due to the dopant ions, inho- 
mogeneous pairing interactions etc. Superconductivity 
arising out of inhomogeneous pairing interactio ns are 
also found in many other systems. lAndersoiI ^ sug- 
gested the possibility of negative-?/ centres in semicon- 
ductors. There are reports of existence of superconduc- 
tivity in silicon based nanostructures with negative-C/ 



centres Negative-?/ models have been used to describe 
the physics of doped bismuthatesi^ A material of more 
recent interest, Tl-doped PbTe, is believed to have a dis- 
tribution of negative-C/ centers i^>^ 

As indicated briefly above, there are two factors that 
lead to an inhomogeneous superconducting state. The 
first one is one body disorder; when one body dis- 
order is large, it tends to localize electrons, and in 
this sense "competes" with superconductivityi^ Effect 
of one body disorder on superconductivity has been ex- 
tensively studied using models and methods of different 
sophisticationi^^i2^i^ The second factor that contributes 
to inhomogeneous superconductivity arises in situations 
where the pairing interaction (such as negative-C/ cen- 
ters) responsible for superconductivity is itself spatially 
inhomogeneous. Systems with a distribution of negative- 
U centers are known to give rise to the "charge Kondo 
effect"—; the material Tl-doped PbTe is believed to be 
one suchi^ Models with inhomogeneous pairing inter- 
actions have been investigated before 
The model usually studied is the attractive Hubbard 
model with inhomogeneous interactions (AHII) is de- 
scribed by the Hamiltonian: 



H = -t 



E 

<lj>cr 



.(1) 



where i,j are site indices of a lattice, t is the hopping am- 
plitude, a is the spin index, cl^ is the electron operator 
which creates an electron of spin a at site i, rii^ = c\^Cic 
is the number operator at site i of spin cr, ^ is the chem- 
ical potential. The interaction C/^ > is site dependent; 
a fraction p (here called the active fraction) of sites have 
C/i = C/ > 0, while C/; = for the other fraction (1 — p) 
of sites. These sites with Ui ^ can be arranged period- 
ically or randomly. It is known that for a given C/, there 
is a critical value Pc of p below which superconductivity 
is kiUed^^i^ 

More recently the above model ([T]), motivated by the 
electronic inhomogeneities in correlated materials, has 
been subjected to extensive numerical simulations i^^i^i^^ 
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FIG. 1: (color online) T = phase diagram in the p — U 
plane. Plots show contours of constant pairing $. The pairing 
amplitude vanishes on the "left" side of the curve marked as 
Pc- (a) Bethe lattice with filling n = 0.44 (b) Flat band with 
filling n = 0.125. 




FIG. 2: (color online) T = phase diagram in the p — fi 
plane. Plots show contours of constant pairing $. The pairing 
amplitude vanishes below the curve marked as pc- (a) Bethe 
lattice with U = 0.6 (b) Flat band with U = 1.0. 



A Bogoliubov-dc Gcnncs mean-field (BdGMF) approach 
is used to obtain the ground state and finite temperature 
properties; these calculations involve averaging over sev- 
eral different [/-disorder realizations. Quantum Monte 
Carlo calculations^! have also been performed on a two- 
dimensional square lattice. Such calculations are numer- 
ically intensive, and attention has been focussed on par- 
ticular values of electron fillings and interaction param- 
eters [/, and the active fraction p. It is useful to have 
a "phase diagram" of the AHII, particularly to compare 
and contrast different experimental systems, and to ob- 
tain regions in the parameter space where interesting 
physics may be expected. Calculation of the phase di- 
agram within the BdGMF approach can be quite time 
consuming; it is therefore desirable to generate the phase 
diagram by means of a simple approach to understand 
its overall structure. 

Motivated by the above discussion, the phase diagram 
of the AHII model is obtained in this paper using a dy- 
namical mean field theory like approachi^ Effects of in- 



teraction strength, active fraction and electron filling are 
investigated systematically. It is found that there are re- 
gions of the phase diagram where "unconventional" su- 
perconductivity is supported. 

The paper is organized as follows. In the next section 
(Section [n|, we discuss the single site formulation of the 
problem. Results of the calculations are presented in Sec- 
tion |TTT1 In the final section (Section lIVp the results are 
discussed, and aspects of the phase diagram not captured 
by the present treatment, including a more speculative 
phase diagram, are discussed. 

II. FORMULATION 

The AHII Hamiltonian given in ([T]) is treated here 
within a dynamical mean field theory like approach. 
A dynamical mean field theory treatment of Hubbard 
like Hamiltonians with one body disorder is treated by 
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iJanis and Vollhardt^ : the present work deals with the 
case where the interaction term is disordered. A single 
site, that hybridizes with an electron bath^ described by 
a bath Green's function Q (this is a matrix in the present 
formulation, see below), represents a typical site of the 
lattice. The imaginary time action for this site is written 
as 

Jo Jo 

- r dTVr^{rWi{r)i^l{rH'^{r) (2) 
Jo 

where [3 = 1/T (T is the temperature), ipa are the Grass- 
mann variables of the site electrons, 

are the Nambu matrices, Q is the matrix Green's function 
(which incorporates the chemical potential fi), V is the 
attractive interaction. The Grassmann variables satisfy 
the Fermionic condition ipa{f3) = — V'(t(0). The interac- 
tion y is a random variable which here is distributed 
according to the probability distribution 

P{V)^{l-p)6{V)+p6{V-U) (3) 

where p is the active fraction, i5(-) is the Dirac delta func- 
tion, and U is the attractive interaction strength at the 
active sites. The disorder averaged partition function can 
now be written as a disorder averaged path integral 

Z = JdVPiV) yp[C,V^.]e-^[**'*^^l (4) 

The site Green's function (expressed in terms of Matsub- 
ara frequencies itOn) is obtained as 

G{iuJn) = (1 -p)giiuJn) +PiG^^{i^n) - S(7 («t^,i ) ) (5) 



where Tjjj is the self energy obtained from the solution 
of the quantum impurity problem with the bath Green's 
function Q and an attractive Hubbard interaction U at 
the site. The Green's function (O represents the lattice 
Green's function as seen from the single site formulation. 
The site self energy is now given by 



'S{iuJn)=G ^{iuJn)-G ^{ioJn)- (6) 



Using the dynamical mean field theory ansata^ that 
the the self energy is "momentum (energy) independent" , 
we obtain the lattice Green's function as 



^{iLUn)= / deg{e) {iLUnl-^T;:-'E{iuJn)) ^ (7) 



where 1 is a 2 x 2 unit matrix, is the Pauli z-matrix, 
£^ ^ e — fi, and g{e) is the bare density of states of the 
lattice. The dynamical mean field theory self consistency 
condition is now obtained by insisting that 



g-\iu;n) = ^^HiUJn) + S(iw„). (8) 
This condition is obtained by demanding that the site 
Green's function as calculated from the quantum impu- 
rity formulation ^ is same as the lattice Green's func- 
tion calculated via ([7]). 

The following paragraphs contain a description of the 
approximate solution, based on the saddle point method, 
of the quantum impurity problem that is used in this 
work. A Hubbard-Stratanovich field A(r) is introduced 
to decouple the interaction term in the particle-particle 
channel. The partition function becomes 



Z = JdVP{V) J V[A*,A]e-^o <^r\Airr J ^^^.^^^^^-S[^\^,V,A] (9) 

where 

5[**,^',V^,A] -y j drdr'**(T) (^g-i(T-r') + \/7A(T-T')j *(t') (10) 



with A defined as 

A( 

The Fermionic path integral is easily evaluated, and the partition function becomes 



jdVP{V) Jv[A*,A]e~i^o d^\Mr)i'~^-<i'^i~is-'+VvA)])^ ^^2) 
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The averaging over the probabiHty distribution can be performed exactly to obtain 

Z ^ J X>[A*, A]e-^^ 

where 



dT|A(7 



-In 



(i-p) dcti-g- 



pdet(-(g- 



UA)) 



(13) 



(14) 



The partition function is now evaluated by introduction of the saddle point approximation which amounts to treating 
A as independent of the imaginary time^; in the present context this approximation is equivalent to the BCS mean- 
field decoupling of the interaction term. With the assumption that A is real (equivalent to picking a particular phase 
of the resulting superconductor), the value of A is obtained by minimizing Sa leading to the equation 



A = 



pe 



lndet[-(e-^+%/[7A)] 



(1 - p)elndct[-e-i] _^pglndct[-(g-i+/l7A)] 

where ^ {Q-^ + VUA)-^). 



XI (^12 (*^") + Gilii'^n)) 



(15) 



Within this approximation the self energy "Eu in ([5|) is 
equal to VU A, where the A obtained from the solution 
of HSl) is used in (fTTj) . 

In present formulation within a dynamical mean field 
theory framework, the saddle point approximation is the 
simplest possible "impurity solver" . The formulation 
based on the Hubbard-Stratanovich fields is amenable to 
more sophisticated, and obviously more computationally 
intensive, treatments such as the Hrisch-Fye quantum 
montc carlo method.— The saddle point approximation 
for the impurity solver is similar to the coherent potential 
approximation!^ 

In the framework developed here, ^/U A has the natu- 
ral interpretation of the disorder averaged pairing gap. 
It should also be noted that there is a possibility of 
introducing a second Hubbard-Stratanovich field in the 
particle-hole channel, which in effect is equivalent to in- 
troducing an additional Hartree potential in the saddle 
point approximation. This extra Hartree potential can 
now be absorbed into the definition of the chemical po- 
tential ^. 

For a given value of U and p, the value of A is cal- 
culated as follows. A typical calculation starts with an 
assumed value of A, and a new value of A is calculated 
using The site Green's function ([5]) and the self 

energy ([6]) are calculated using the new value of A. The 
site self energy is used in ([7]) to obtain the lattice Green's 
function, and a new bath Green's function is generated 
using ([5]). This process is carried out until the values of 
A, self energy S are within a specified tolerance of each 
other in two successive iterations. All self-consistency 
calculations are done at fixed chemical potential /i and 
the number of electrons n is obtained after convergence 
is obtained. The chemical potential is then adjusted to 
that the number electrons is obtained to be the desired 
value. The calculations reported here are performed by 
evaluating all Matsubara sums as integrals along the real 



frequency axis, and the self consistency condition also 
enforced on the real frequency axis. 

III. RESULTS 

In the single site formulation presented in the last sec- 
tion, the information regarding the lattice enters the for- 
mulation only via the density of states 5(e). Since the 
goal of this paper is to understand the generic features 
of the phase diagram of the AHII model, densities of 
states with simple analytical forms that capture some 
features of the real lattice systems are adopted. Two 
cases are considered. The semicircular density of states 
corresponding to a Bethe lattice^ with 

2 

g{e) = -Vl -e^ -1 < e < 1 (16) 

TT 

and the flat band density of states 

5(e) = -l<e<l. (17) 

Energy is measured in the units of half bandwidth of 
the systems, and hence the condition — 1 < e < 1 in 
both cases. Other parameter values of quantities such as 
U, fi, T are all henceforth dimensionlcss ratios of these 
quantities and the half bandwidth. 

Superconductivity is monitored by computing the pair- 
ing amplitude <&: 

$ = (i/^fV^) = -- / dc^3Gi2(c^) (18) 
^ J —00 

Vanishing of $ implies absence of pairing and supercon- 
ductivity. Clearly, existence of a nonzero value of $ auto- 
matically does not imply global superconductivity. This 
point is discussed in more detail later in the paper when 
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FIG. 3: (color online) Spectral function for the Bethe lattice 
with U — 0.8, n = 0.44. The four panels show the spectral 
function for increasing values of the active fraction p. 

the Bose-Einstein condensation like phenomenon in such 
systems is discussed. 

Results at zero temperature are presented first followed 
by results at T > 0. 

A. Results: T = 

For a given density of states, the phase diagram is de- 
termined by three parameters: the interaction strength 
U, the active fraction p and the filling n. Fig.[T]shows the 
phase diagram of the system in the p-U plane for different 
fillings. For a given filling n we see that there is a range of 
interaction strength U for which there is a critical value 
of the active fraction that is required to produce a 
nonzero pairing amplitude In the case of the Bethe 



lattice the value of Pc decreases with increasing filling n, 
while for the flat band case, Pc is essentially insensitive 
to filling. This result is a reflection of the fact that pc 
is affected by the bare density of states (at the chemi- 
cal potential); the bare density of states increases with 
increase in filling (up to n < 0.5) for the Bethe lattice, 
and hence the decrease of pc- The critical active fraction 
Pc is insensitive to filling in the flat band case since the 
density of states is constant. This is more clearly illus- 
trated in Fig. [2] which shows the phase diagram in the 
p — n plane for two values of U. It is evident that Pc in- 
creases with decreasing n, and this increase is related to 
the decrease in the bare density of states (see Fig. [2ja)). 
Interestingly, the sensitivity of Pc on n decreases with in- 
creasing interaction strength U. The present calculation 
reveals an interesting new result. For values of U larger 
than a critical value Uch (this value depends on the fill- 
ing, i.e., bare density of state at the chemical potential), 
the critical value of the active fraction pc becomes van- 
ishingly small (see Fig.[T|). Thus if C/ > Uch, even a small 
concentration of impurities can produce a non-zero pair- 
ing amplitude. In the same vein, there is another critical 
value of the interaction Ucl- If the interaction strength is 
below UcL (which, again, depends on the bare density of 
states at the chemical potential) even a small dilution of 
the active fraction from unity kills the pairing amplitude! 
The spectral function, 

iV(cj) = --5Gii(w), (19) 

TT 

also provides interesting information regarding the na- 
ture of the electronic state. Fig. [3] shows spectral func- 
tions for the Bethe lattice with U = 0.8, n = 0.44 (re- 
sults for the case of the flat band qualitatively similar), 
for various values of the active fraction p. The spectral 
function at p = 1 has a gap with a characteristic BCS 
(Bardeen-Cooper-Schrieffer) singularity^ in the spectral 
function near the gap edges. On the other hand for p < 1 
it is seen that the singularity is "smeared out" , by ap- 
pearance of "mid-gap states". The calculation suggests 
a possible spatial distribution of gaps in the system with 
different regions of the lattice developing different gaps. 
It is, of course, not possible within the present frame- 
work to study the gap distribution, but further detailed 
simulations could throw more light on the nature of the 
inhomogeneous state. 

B. Results at T / 

Fig. [4] shows the phase diagram of the AHII model 
in the p — T plane. Several interesting features can be 
seen. For all active fractions with a non-vanishing pair- 
ing amplitude at zero temperature, there is a tempera- 
ture Tc at which the pairing amplitude (and A) vanishes. 
There are three regimes of active fraction that give rise 
to very different flnite temperature phenomenon. For ac- 
tive fractions just above the critical value (marked pc in 
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FIG. 4: (color online) Phase diagram in the p — T plane. Plots show contours of constant pairing $. The nearly straight line 
marked Tc represents a continuous transition. The attractive interaction U = 0.8 in all cases, (a) Bethe lattice with n = 0.072 
(b) Bethe lattice with n = 0.44 (c) Flat band with n = 0.125 (d) Flat band with n = 0.45. 



the panels in Fig. [4]), the transition to a regime of van- 
ishing pairing amplitude takes place by an abrupt (first 
order) transition (see, for example, region 0.2 < p < 0.6 
in Fig.[4ja)). On increase of the active fraction, a second 
regime appears (see, for example, region 0.6 < p < 0.75 
in Fig.H]) where there are two transitions. In this regime, 
with increase of temperature from zero, there is a first 
order transition where $ (and A) undergoes a sudden 
jump and obtains a smaller non zero value. With further 
increase of temperature the pairing amplitude vanishes 
continuously to zero. Interestingly, the first order line in 
the p — T plane appears to end at a "critical point" (such 
as that marked by C in Fig. Ufa)). With further increase 
of the active fraction, a third regime is attained (j> > 0.75 
in Fig. [4{a), for example), where Tc depends essentially 
linearly on p, and the transition is continuous. These ob- 
servations can be clearly seen by a study of Fig. [5] which 
shows a plot A(T)/A(0) as a function of T/Tc, where the 
three types of behaviour arc shown. Indeed, the features 
are generic and do not appear to depend on the shape of 
the bare density of states; they are clearly seen in both 



the Bethe lattice and flat band cases. The "sizes" of the 
three regimes are, however, strongly affected by the elec- 
tron filling n. This is most clearly seen in the fiat band 
case where the critical value of pc is insensitive to elec- 
tron filling. However, the finite temperature behaviour 
strongly depends on electron filling - compare Fig. [D^c) 
and Fig. IDJd)), in the latter case the second regime of 
active fraction with two finite temperature transitions is 
strongly suppressed. 



IV. DISCUSSION AND CONCLUSION 

This section contains a summary of the results ob- 
tained in this paper, and a discussion of the full phase 
diagram of the AHII model. The present calculation of 
the phase diagram is based on a single site dynamical 
mean field theory like approach. The calculation shows 
for a certain range of the strength of the interaction pa- 
rameter U which depends on the bare band structure 
and electron filling, there is a critical active fraction pc 
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FIG. 5: (color online) Variation of A with temperature T for 
the flat band U — 0.7, n = 0.45. The four panels show the 
behaviour of A for increasing values of the active fraction p. 
The critical value of the active fraction for this case pc = 0.62. 



below which there is no pairing amplitude. However, for 
large enough values of the interaction parameter, even 
an infinitesimal value of the active fraction is sufficient 
to produce a non vanishing pairing amplitude. For active 
fraction p greater than pc the electron spectral function 
shows a features with "mid-gap states", and the BCS 
singularity of the spectral function is smeared out. 

For cases which show a nonvanishing pairing ampli- 
tude $, three types of finite temperature behaviour arc 
found. For p close to Pc, there is a discontinuous transi- 
tion at a finite temperature, and for p close to unity, there 
is continuous transition (BCS like behaviour) to a state 
without pairing amplitude. There is a intermediate range 
of active fractions, where the transition to a non paired 
state takes place in two steps - "non-BCS behaviour". 
As the temperature is increased, there is a first order 
transition to a state with smaller Further increase 
of temperature causes a continuous transition to a state 
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FIG. 6: (color online) Schematic phase diagram of AHII. 
The region indicated by "Non-BCS superconductor" indicates 
the unconventional superconducting behaviour found in the 
calculations presented in this paper. The region indicated 
by "Heavy Fermi Liquid " can be inferred from the work of 
iTaraphder and Colemar^ . The region with "large t/" and 
large p resul t in a "BEC of pairs" as is known from the work 
of Ra nderia^ . and from the estimate of superfluid stiffness 
given in fig. [T] The region with large U and small p is specu- 
lated to have a "pair glass" ground state. 



with no pairing amplitude. It is tempting to speculate 
that the state attained up on the first order transition 
has a pairing amplitude, but no superconductivity. The 
physical picture of such a state is that of "puddles of elec- 
trons" with non-zero pairing amplitude without a global 
phase necessary for superconductivity. Such a state is 
likely to show "psuedo-gap" like features, for example, a 
reduced spin susceptibility. Clearly, this finding of the 
present calculation needs more attention, and the region 
of the phase diagram where this phenomenon is found 
needs further detailed investigation. It is interesting to 
note that calculations based on BdGMF— also show a 
regime of U and p which show anomalous behaviour of 
A as a function of T . 

The present formulation is based on a single site for- 
mulation and averages over all the spatial correlations. 
However, as noted above, a very interesting region in the 
phase diagram is revealed, and suggests possibility for 
further investigation. Further, the approximate treat- 
ment based on the saddle point approximation does not 
include quantum fluctuations. It is believed that the in- 
clusion of these quantum fluctuation effects are not likely 
to change the qualitative features of the present single 
site calculation; this is suggested by the iterated pertur- 
bation theory based dynamical mean flcld theory of the 
attractive Hubbard model4^ It must be noted that most 
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FIG. 7; (color online) Dependence of the superfluid density 
ps on the interaction strength U in the Bethe lattice with 
n = 0.44 and p = 0.4. 

of the previous work cited above arc based on two dimen- 
sional systems, mostly square lattices. Long wavelength 
fluctuations, crucial in two dimensional systems, cannot 
be treated within the present framework. The present 
work, therefore, is more applicable to higher dimensional 
systems such as negative- ?7 centre systems etc. 

The paper is concluded with a discussion of the com- 
plete phase diagram of the AHII model. Based on 
the calculations presented here and on published results 
quoted earlier, the nature of the electronic state in dif- 
ferent regions of the parameter states can be inferred. 
In the regions with nonzero pairing amplitude (as ob- 
tained from the present calculation), different types of 
electronic states may be found as indicated in Fig. [HI For 
intermediate U (compared to the bare bandwidth) and 
large p a BCS superconductor is obtained. On the other 
hand for smaller value of the active faction p, a "non- 



BCS" superconductor is obtained (as discussed above), 
with the possibility of a high temperature "pseudo-gap" 
phase. For larger values of the interaction strength there 
is a crossover from BCS to BEC like behaviour where the 
electrons form pairs and Bose condense4^ In the present 
calculation, this behaviour is inferred by the calculation 
of an estimate (based on the kinetic energy) of the super- 
fluid densityi^ As shown in Fig. [71 the superfluid density 
Ps falls with increasing U indicating a crossover from BCS 
to BEC (Bose-Einstein condensation) behaviour. The 
effect of the random attractive interaction in this BEC 
regime needs a more careful investigation than that given 
here. For low values of the active fraction and interme- 
diate values of the interaction strength, the most likely 
ground state is a heavy Fermi liquid^! engendered by the 
charge Kondo effect. It is also possible that for small fill- 
ings, large interaction strengths, one could obtain a pair 
glass, where electrons are localized at the negative U cen- 
ters. Clearly, the "boundaries" of the regions indicated 
in the phase diagram will be determined by the third fac- 
tor in the problem, namely electron filling. The nature of 
the electronic states in different parameter regimes of the 
Ami model does resemble various systems discussed in 
the introductory section. It will be interesting also to ex- 
plore the possibility of a direct experimental realization 
of the AHII model in cold atom optical lattices ^li^i^ 
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